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Definition 

(i)th 4x4 Link Transform of Mechanism D 
Mechanism Jacobian matrix 
Positive Definite diagonal scaling matrix 
Transpose of matrix M 
Inverse of matrix M 

Product Aj)j for all i in the (j)th time interval 
Gravity gradient vector 
Four dimensional position vector 
Scalar or dot product of two vectors 
Vector or cross product of two vectors 
Norm or length of a vector 
Parallel condition between vectors a and b 
Vector composed of mechanism joint angles 

(i) th Mechanism joint angle 
Ideal mechanism joint solution 

(j) th Estimate of ideal mechanism joint solution 


UNUSUAL TERMS 


downhand welding - Welding with the part maintained in an orientation such that the weld 
puddle lies in a horizontal plane. 

GTAW - Gas tungsten arc welding is performed with a tungsten electrode in an inert gas 
atmosphere. This process is also call TIG or tungsten inert gas welding. 

off-line programming - A method of specifying the motion path of a robot in an external 
computer. 

wire feed - When welding with a consummable electrode, the electrode is usually supplied 
to the welding process in the form of a wire fed into the arc directly in front of the 
electrode as it moves along the weld seam. 


A GENERALIZED METHOD FOR AUTOMATIC DOWNHAND AND WIREFEED 
CONTROL OF A WELDING ROBOT AND POSITIONER 


INTRODUCTION 

Automated arc welding imposes several requirements on robot/positioner controls 
that are unique to many other applications. The most important of these are: (a) the 
desirability to weld in such a position that the effect of gravity on the molten weld pool is to 
aid the metal fusion process rather than to impede it, (b) the need to maintain constant 
travel speed of the torch relative to the part to insure a constant heat input to the weld, (c) 
for Gas Tungsten Arc Welding (GTAW) processes, the need to maintain proper weld wire 
orientation, and (d) the need to accurately track the weld joint to insure proper fusion at all 
points along the weld. The latter requirement stems from the need to adapt to variations 
between the joint and the pre-programmed path of the manipulator due to thermal 
distortion, dimensional variations in the position of the weld joint on the welded assembly, 
dimensional variations in the weld assembly, variations in the weld fixturing, and variations 
in the equipment (i.e., the manipulator and the positioner). Cook [1], Richardson, Farson, 
Jones, and Rogers [2], and others have addressed the joint tracking problem from the 
standpoint of both optical and electrical arc sensing. 

The purpose of the algorithm described in this paper is to address the other three 
requirements of an automated arc welding system, i.e., welding in the downhand position, 
maintaining weld travel speed, and in GTAW proper weld wire orientation. A 
transformation algorithm is derived that permits programming the weld path teach points 
with the assembly in any convenient position. It is only necessary to position the torch at 
each programmed position such that it is at the desired orientation relative to the weld 
joint. The transformation algorithm computes the required kinematic transformation 


matrix for the manipulator, the manipulator velocity, the position vector for the positioner 
axes, and the velocities of the positioner axes to insure that the weld is performed in the 
downhand position and that the weld travel speed is constant with wire feed orientation 
controlled. The algorithm is applicable to any arc welding robot/positioner system. It is 
only necessary to be able to download the kinematic transformation matrix, positioner axes 
values, and velocity information to the robotic controller. 

The weld process requirements are discussed first. This is followed by a derivation 
of the algorithms, results of a numerical simulation coupled to the graphics simulation, and 
a discussion of an implementation of the algorithms in hardware. Applications are 
described for welding of the Space Shuttle Main Engine. 

DOWNHAND WELDING 

For most arc welding processes the highest deposition rate can be obtained when 
welding is performed in the flat (or downhand) position. The downhand position is defined 
as that position in which welding is performed from the upper side of the joint and the face 
of the weld is approximately horizontal [3]. The highest deposition rate is achievable in the 
flat position because gravity acts to keep the molten metal in the joint. This permits the use 
of higher welding current, with a resulting larger molten weld pool. Studies of flat position 
welding versus welding in other positions have proven that a heavy weldment can be made 
in approximately one-half the time if all of the welds are made in the flat position [4]. 

Another reason for positioning the workpiece for downhand welding is to improve 
weld quality. The higher deposition rates achievable with downhand welding mean fewer 
passes are required to complete a multi-pass weld. This in turn can be expected to reduce 
the amount of distortion that may be encountered. Also, the desired weld bead contour is 
much more easily maintained in the flat position than in other positions. For example, a 
fillet weld is defined as a weld of approximately triangular cross section joining two 
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surfaces approximately at right angles to each other [5]. In the downhand position, each 
surface would be at an angle of approximately 45-degrees with respect to the horizontal. 
The face of the weld, corresponding to the hypotenuse of the triangular cross section, 
would be horizontal. In the flat position, gravity would act to maintain this desired weld 
bead profile. On the other hand, if one of the surfaces being welded lies in the horizontal 
plane then the other surface would be in the vertical plane and the face of the weld would 
be undesirably at an angle of 45-degrees with respect to the horizontal. This position of 
welding is termed the horizontal position [5], The effect of gravity is now to cause the 
molten weld pool to sag toward the horizontal surface. This may result in an unacceptable 
weld profile and possible under-cutting in the vertical surface. 

POSITIONING OF WELD ASSEMBLY 

Two types of positioner controls are typically used for automatic, robotic arc 
welding. With one type, the weld assembly is indexed to a programmed weld position, and 
the weld is performed by moving the manipulator only. When the weld is completed the 
assembly is indexed to the next programmed weld position, and the next weld is performed. 
This scenario is repeated for all welds (of relatively straight contour) on the assembly. For 
circular welds the assembly is indexed to a programmed welding position and then rotated 
while the weld is performed. Even though time is wasted during indexing intervals (when 
welding is not performed), this type positioner control is adequate for assemblies 
containing relatively short, straight weld contours or circular contours. The programmed 
welding positions may be chosen such that the welding is performed in the downhand 
position. Programming is relatively easy with a teach pendant. Constant weld travel speed 
(and hence constant heat input) is easily maintained because the weld contours are 
assumed to be straight or circular and because the positioner and the manipulator are not 
in motion simultaneously. 
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With the other type of positioner control the motion of the positioner is coordinated 
with that of the manipulator such that both can be in motion simultaneously. This type of 
control is particularly important for assemblies containing irregular joint contours. 
Typically, the control algorithms (used with commercially available robots) permit 
manipulator motion from one programmed point p(i) to a second programmed point 
p(i+l) along a straight line (relative to robot base coordinates) at a constant velocity 
[g(i+ l)-p(i)]/T. The velocities of the positioner axes are automatically set to be completed 
in the same time interval T. This capability within itself does not guarantee either 
downhand welding or constant weld travel speed. However, the former can be obtained by 
using the teach pendant to position the weld assembly (for each point to be programmed) 
such that the weld is in the flat position at that point. At each such point the manipulator is 
located (via the teach pendant) to provide correct position and orientation of the weld 
torch. This is a tedious task and does not ensure that the welding travel speed will be 
constant. Indeed, the welding travel speed will not be constant on any but the simplest of 
weld contours. This is due to the relative motion (between the assembly and the torch) that 
is not taken into account in establishing the manipulator travel speed. GTA welding further 
complicates the robot programming problem since filler wire must be fed into the arc 
plasma in a controlled manner. For most applications the filler wire is fed directly in front 
of the moving plasma. As stated in the introduction, these problems (i.e., maintenance of 
downhand position, constant travel speed measured with respect to the weld joint and 
proper weld wire orientation) are overcome with the transformation algorithm discussed in 
this paper. 


OVERVIEW OF ALGORITHM DEVELOPMENT 
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The algorithms described use an iterative method to determine the robot and 
positioner joint angles needed to traverse a weld path program, maintain the orientation of 
the part so that welding in the downhand position, and maintain weld wire orientation. This 


functionality is achieved in three steps by the development of three algorithms each of 
which achieved more control of the robot/positioner than its predecessors. The first 
algorithm maintains torch to part position/velocity control using the eight degrees of 
freedom of the mechanism formed by combining the positioner and robot. The additional 
constraint of minimizing the weighted sum-of-squares of joint angle displacement is used to 
eliminate the redundancy inherent in the combined eight DOF mechanism. The second 
algorithm achieves torch position/velocity while simultaneously controlling the orientation 
of the part so that welding is performed in the downhand position. The second algorithm 
differs from those previously published [6,7,8] in that now robot and positioner are treated 
as a single mechanism. The third algorithm is similar to the second with the added 
constraint of maintaining weld wire feed orientation. 

DEFINITIONS 

We begin our description of the algorithm by stating several definitions. The 
homogeneous transformation is commonly used in the robotics literature [9,10] and consists 
of a 4 by 4 matrix, J'^Aj or more simply Aj, as shown in equations (1) and (2). 
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The upper left 3 by 3 partition of Aj is a rotational transformation of a vector space, 
j, whose orthogonal basis are the vectors n, p, and a, with respect to the j-1 frame. Vector p 
defines the displacement of the origin of frame j with respect to frame j-1. In defining the 
motion of a robot reference, frames are defined fixed within each link [11] and the position 
and orientation of an N link robot’s hand reference frame is defined in the base coordinate 
frame by Tr in equation (3). 
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[T R (i)] 


(3) 


A R1 A R2 A R3* * * A RN 

Each link transformation, Aj, is a function of the (j)th joint angle for a revolute joint 
or the (j)th extension for a prismatic (sliding) joint. A similar transformation may be 
defined for the robot’s M DOF positioner given in equation (4). 

[Tp(i) ] = A Pl A p2 A P3 * * * A PM (4) 

ALGORITHM I: TORCH POSITION CONTROL 

The position of a tool center point (in our case a welding torch) may be defined by a 
series of transformations which includes Tr. Figure 1 depicts a transform graph of our 
robot workcell. The origin of our workcell is designated by point O in the figure. Zr, a 
fixed transformation, defines the position of the robot’s base coordinate frame with respect 
to the workcell frame. [Tp(i)j is the robot transformation as defined in equation (3) with 
the variable, i, indicating the (i)th programmed point in our welding path. Finally, the 
homogeneous fixed transformation, E, is the end-effector or torch transformation. In 
defining E in the same manner as equation (1) the vector, a, is chosen to be along the torch 
plasma flow and the vector, g, is the point in the torch plasma where contact with the seam 
is required. 

In a similar manner the position of the part weld path may be defined by a series of 
transformations: Zp is a fixed homogeneous transformation which translates the part 
positioner’s base reference frame to the workcell frame; [Tp(i)j, the positioner 
transformation matrix defined as in equation (4); [G], the part frame to positioner 
mounting frame transformation; and [P(i)J, the transformation which relates the desired 
position of the weld torch reference frame (at the (i)th program point) to the part 
reference frame. The transform, [P(i)J, may be determined either "on-line" by using a teach 
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pendant to move the robot through the desired torch positions (not necessarily maintaining 
the downhand position), or [P(i)] may also be determined by some "off-line" technique 
(e.g., interaction with a CAD/CAM system). 

The condition that the welding torch is located at the desired path program point, 
represented by the matrix [P(i)], is met when equation (5) is satisfied. 

z R [T R (i)]E = Zp[T p (i)]G[P(i)] (5) 

Re-arranging equation (5) to form equation (6) and referring to figure (1) we note 
that the desired torch position, [P(i)J, is equal to the actual torch position when there is no 
error. 


P(i) = [G- 1 ][Tp(i)" 1 ][Z p ” 1 ][Z r )tT r (i)]E (6) 

Expressing equation (6) in terms of our link transformations we get equation (7). 
(Note: For this example the robot has six DOFs and the positioner has two DOF’s.) 

P(i) = (G- 1 ][A 2p - 1 ]tA lp - 1 ][Z p - 1 ]Z r A lr A 2r A 3r A 4r A 5r A 6r E 

( 7 ) 

An error in torch position may be expressed by defining the transformation, [Err], as 
the transformation between the desired torch position and the actual position. Equations 
(8) and (9) represent this more general condition. 

P(i) = [Act] [Err] (8) 
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[Act] — [G [^ 2 p 1] [^lp [2p ■*"] Z r A^ r A2 r A3 r A4 r A5 r A6 r E 

( 9 ) 

Figure (2) depicts our modified workcell transform graph in which the torch position 
error is represented. This error may be determined easily from the command, [P(i)], and 
the [Act] matrices as shown in equation (10). We may also represent the error, [Err], as a 
translation and three rotations about the x,y and z axes shown in equation (11). 

[Err] = [Act” 1 ] [ P ( i ) ] (10) 

[Err] = [Trans (dx, dy, dz) ] [Rotx(d0 x ) ] [Roty(d0y) ] [Rotz (d0 z ) ] 

( 11 ) 

For very small rotational errors, d<p x , d<f>y and d<p z , the order of rotation becomes 
less distinct. This may be shown by expanding equation (11) in terms of sines and cosines 
and using the small angle approximations: cosine(d0) is equal to one minus d 0; and 
sine(d0) is equal to d<p. Terms involving higher orders of d<p are then eliminated resulting 
in the small angle approximation of [Err] in equation (12). 
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( 12 ) 


Equation (12) may be re-written from the matrix form into a vector form as shown 
in equation (13). A vector, v, which is transformed by the matrix, [Err], is equivalent to the 
same vector affected by a displacement, v^, and rotation about $_^ rx with v^j. and < p C TT 
defined in equations (14). 
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[Err ] v = 


( 13 ) 


[I]v + Vgj» r + 0 err X v 



dx 


d0 x 


dy 


d0 y 

V = 

—err 

dz 

' ^err “ 

d0 z 


1 


0 


An examination of the matrix form in equation (10) and the vector form in equation 
(13) reveals that as the error displacements and rotations vanish [Err] becomes the identity 
matrix and figure (2) becomes equivalent to figure (1). 


A six degree-of-freedom error function, fi>p_ err , may be formed by concatenating 
Vgjj. and Since control of the torch motion will be related to the path, [P(i)], defined 
in the part frame, fi>p_ err is shown in equation (15) transformed to the part reference frame 
by [Act], the homogeneous transform of the actual torch position. Matrix [Q] in equation 
(16) removes the 4th, scalar component, from the normalized homogeneous vector 
reducing it to a 3D vector. The resulting 6D vector is given in equation (17). The 
independent variable is an eight-dimensional vector, 0, in equation (18). The components 
of 0 are defined in terms of the positioner and robot angles in equation (19). This 
definition is consistent with the definition of the robot/positioner combination as a single 
serial kinematic chain defined from the part to the end-effector. 


[Q] [Act] v err 
[QJ [Act 3*^ 


( 15 ) 


where, 


[Q] = 


10 0 0 
0 10 0 
0 0 10 


( 16 ) 
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— p-err = [dx,dy,dz,d0 x ,d0 y/ d0 z ] 


— p-err (— ) 


(17) 


T _ 

0 = [ 0 ^ i ®2 • ®3 • ® 4 • ® 5 ' ® 6 > ® 7 ' ® 8 ^ 

® = [ e 2P' e lP' e iR' e 2R' e 3R' e 4R' 9 5R' e 6R3 T 


(18) 

(19) 


The vector error function has the property that the norm of £6p_ err vanishes as the 
torch position and orientation errors vanish (equation 20). This condition is also expressed 
in equation (21) when the function approaches the vector, 0, a vector having six zero 
components. 


[(f6 p _ err ) T ][f6 p _ err ] -> 0, as [Err] -> [I] (20) 

— p-err (—) _> £ (where 0 = [0, 0, 0, 0, 0, 0] T ) (21) 

Newton’s method has been used by Whitney [13] to solve equation (21) for 0. The 
iterative equation implementing the Newton’s method is shown in equation (22). The 
predicted solution vector, 0j + j, is equal to the previous solution vector, 0j, modified by a 
correction term as shown. One factor in the correction term of equation (22) is the 
Jacobian matrix which is defined in equation (23). (Equation (A-3) of Appendix A defines 
the Jacobian matrix in terms of the robot and positioner transformation matrices.) Since 
the Jacobian matrix, J, has eight columns corresponding to the eight joint angles of our 
combined mechanism and J has only six rows for our differential translations and rotations, 
J is non-square and the pseudo-inverse of J must be used (see equation 24). The scalar 
factor, h, in equation (24) is the iteration step size. The use of the pseudo-inverse results in 
a solution of the under-determined system of equations for the eight component angles of 0 
that satisfy the torch-to-part position requirement and minimize the sum of the squares of 
the eight joint displacements. A derivation of equation (24) is presented in Appendix B. In 
solving for 0, equation (24) is repeated until the weighted vector norm in equation (25) is 
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satisfied. The elements of the diagonal weighting matrix, K, and epsilon may be chosen to 
increase or decrease the tolerance to position and orientation errors in the final solution. 

— j+1 “ — j — p-err^— j)3 ^— ] p-err^— j) 

( 22 ) 

[J] = Cg| f6p-err(&j)3 

( 23 ) 

e j+1 = e j + h* [ j t ( jj t ) _1 3 f6 p _ err (e j ) 

( 24 ) 

tf6p. err (e) T ][K][f6p_ err (0)] < epsilon 

( 25 ) 


ALGORITHM II: POSITION CONTROL WITH DOWNHAND TORCH ALIGNMENT 


In this section we will derive the control needed to maintain the proper weld 
orientation. To accomplish this we will first examine the robot’s path program to determine 
the desired local vertical in the part frame, and secondly we will develop additional 
constraint equations to augment the forcing function defined in equation (15) achieving the 
correct torch-to-part position while maintaining proper part orientation for downhand 
welding. 

The robot’s path was defined above as a series of homogeneous transformations, 
where [P(i)] denotes the desired position of the torch defined in part coordinates for the 
(i)th program step. The first three columns of (P(i)] are the basis vectors of a reference 
frame tied to the torch and the fourth vector is the desired position of the torch frame’s 
origin (see equation 26). 

C P ( i) 3 = [ n(i) o(i) a(i) p(i) ] ( 26 ) 

Figure (3) depicts a segment of the path sequence defined by [P(i)] in the part 
reference frame. In figure (3) a(i) and p(i) are the third and fourth columns of [P(i)]. The 
program points have been chosen so that the secant, g(i), between adjacent points, p(i-l) 
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and p(i), approximates the weld path tangent with sufficient accuracy. A unit vector, n s (i), 
along the vector s(i) is defined in equation (27). In a majority of weld programs the vector, 
a(i), directed along the torch plasma flow, the path secant, s(i), and the desired gravity 
gradient, ng(i), will be co-planar. Additionally, ng(i) will be perpendicular to the path 
tangent which is approximated in this case by the secant vector, s(i). Figure (4) depicts this 
relationship and equation (28) defines the vector ng(i) which is to be aligned with the 
gravity gradient vector in the workcell. n s (i) is the unit vector along §(i). 

E(i) - E(i-l) 

n s d) = -7 : : i“ (27) 

|| E (i) - E(i-l) II 

a(i)-[n s (i) ’a(i) ]n s (i) 

ng(i) = (28) 

II a(i)-[n s (i) *a(i) ]n s (i) || 

Downhand position control is achieved when the part is positioned so that Hg(i) and 
the Earth’s gravity gradient, g, are aligned at the (i)th program point. A rotation vector 
(equation 29) equal to the cross product of ng and g is used as a rotational forcing function 
that will tend to align ng with the local gravity gradient, g. We note that as in the case of the 
positional forcing function defined in equation (15), this new gravity alignment function 
vanishes as the error in alignment diminishes (equation 30). 

fg = Eg(i) X g (29) 

II fg II -> 0 as Eg -> Eg I I g (30) 

The coordinate frame for our robot’s workcell in this study was chosen with the x 
axis vertical, resulting in the gravity gradient, g, being equal to (-1 0 0 0)^. Since g is fixed 
with respect to workcell coordinates the vector forcing function defined in equation (29) 
will lie in the y-z or horizontal plane, therefore the vertical component of the vector forcing 
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function may be eliminated reducing the forcing function to a 2D vector as shown in 
equation (31) with [R] defined in equation (32). 


f 2 — 


err 


where, 


[R] = 


[R] ( ( [Z p ] [A lp ] [A 2p ] [G]ng(i) ) X g) 


0 10 0 
0 0 10 


(31) 


(32) 


The 2D vector error function, £2g_ err , also has the property that its norm vanishes as 
Hg(i) and g are aligned as shown in equation (33). 


[ (f2g_ err ) T ] [£2g_ err ] -> 0 as ng(i) X g -> 0 (33) 

A combined error (equation 34) function, f8 t _ err , is formed by augmenting the two 
previously defined error functions. This function approaches the zero vector, 0, as all 
constraints are met (equation 35). 


— t-err 


— p-err 


f 2 — 


err 


-t-err -> £ 


(where 0= [0,0,0,0,0,0,0,0] T 


(34) 


(35) 


Once again Newton’s method is used to iteratively solve equation (35) for the eight 
joint angles required to control torch position and weld seam orientation. This algorithm is 
given by equation (36). The evaluation of this equation is repeated until the error 
measurement in equation (37) falls below the acceptable tolerance, epsilon. Elements in 
the diagonal matrix [K] are chosen to vary the relative tolerance on torch to weld seam 
positioning accuracy and downhand orientation of the part. If these diagonal elements are 
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denoted by k^, then: torch-to-seam position tolerance is affected by for i = 1,2,3; torch- 
to-seam orientation tolerance is specified in k terms 4,5 and 6; while terms 7 and 8 control 
downhand part orientation tolerance. 

fli+1 = fli + h* [ J -1 ] f8. t _ err (36) 

[(f8 t _ err ) T )[K][f8 t _ err ] < epsilon (37) 

The derivation of the Jacobian matrix, [J], which is needed for the implementation 
of the Newton’s method specified in equation (36) is found in Appendix A (equation A-23) 
while the inverse of this Jacobian is given in Appendix C. 

ALGORITHM III: WIRE FEED 

This section develops the control needed to adjust the wire feed mechanism so that 
wire is fed in front of the moving torch. In our work to this point we specified the torch’s 
position and orientation with-respect-to the part. This specification established 
displacements in x, y, z and orientations of pitch yaw and roll requiring six degrees-of- 
freedom (DOF). The additional constraint equations developed to maintain the downhand 
orientation required two more DOF’s. Since our total mechanism included the two DOF 
positioner and six DOF robot, our system is totally determined. If further constraints on 
torch motion are needed we must either add an additional DOF’s e.g., an auxiliary joint, or 
we may relax a previously imposed constraint. The later approach is used here. Positioning 
of the torch within the weld seam requires that control be maintained in only five 
dimensions: x, y, z, pitch and yaw about the torch tip. The sixth DOF, roll about the torch 
axis, does not affect the welding process and will be used to control wire feed orientation. 

Wire feed control is necessary to maintain the position of the weld wire’s entry into 
the plasma from the front of the torch as it moves along the seam. The method used is 
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similar to that of the previous section. We will define a constraint equation that will force 
our solution to include proper wire feed orientation. To begin our definition of this forcing 
function, fl w _ err we must expand the definitions of the path program vectors presented in 
figure (3). Referring to figure (5) we observe a close-up detail of the path. As in the 
previous definitions, vectors $(i) and §(i+ 1), are path secant vectors with the corresponding 
unit vectors n s (i) and n s (i+ 1), respectively. The vector ng(i) is the desired vertical vector 
used to correct for downhand orientation at the (i)th program point. Also shown in figure 
(5) is n w (i), a unit vector along the direction of the wirefeed and pointing counter to the 
wire feed direction. This vector is fixed with respect to the torch and the end-effector 
reference frame, but it is designated n w (i) in figure (5) since it varies within the part frame 
and is associated with the (i)th path-program point. The end-effector’s approach vector, 
a(i), is shown co-linear to the axis of rotation of the robot’s last joint, g^r. This 
relationship is necessary to insure that a rotation of G^r will affect wire-feed orientation 
without displacing the point at which welding occurs. 

The algorithm proceeds by first solving the inverse kinematic equations iteratively 
using equations (36) and (37) resulting in an arm position that achieves welding in the 
downhand position. The robot’s last joint angle is then modified to achieve proper wire 
feed orientation. Referring to figure (5), proper wire feed orientation occurs when n w (i), 
n s (i) and a(i) are co-planar and the angle formed by the positive directed segments of n w (i) 
and a(i) includes n s (i). 

A scalar forcing function, fl w . err , defined as in equation (38) vanishes when the 
optimum wirefeed position is achieved. The sign of this function is such that it may be used 
iteratively to drive G^r to the optimum position when the magnitude of G^r is defined by 
the right-hand-rule about a(i). 

^w~err — ^ Ilg(i) 1 §(i) (38) 
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This function may be re-written in terms of the kinematic matrices as shown in 
equation (39). 

fl w-err = [ ( CG _1 ] [T p (i) _1 ] [Zp" 1 ] [Z R ] [T R (i) ]En w ) xn s (i)]‘a(i) 

(39) 

In evaluating equation (39) the transformation angles which include the two 
positioner angles, 9jp and 02 p, and the first five robot angles, 0 ^ to 05 ^, are fixed at the 
result derived from the solution of equation (36) and (37) in the previous section on 
downhand control. The robot’s last degree-of-freedom, 0 ^p, is derived by repetitive 
evaluation of equation (40) until the condition in equation (41) is satisfied. The constant, h, 
in equation (40) adjusts the iteration step size. 

e 6R (j+1) = e 6R (j) + h * fl w-err ( 40 ) 

k * (fl w _ err ) 2 < epsilon (41) 


SIMULTANEOUS TOTAL SOLUTION 


An alternative to the separate solution of wire feed orientation and downhand path 
control is a method using a concurrent form of the iterative equations. The concurrent 
solution may be written by combining the equations (36), (37), (40) and (41) into a single 
pair of iterative and constraint equations. A generalized nine dimensional vector forcing 
function, £2 t _ err is defined in equation (42) as the concatenation of our £6p_ err , f2g. err and 
flw-err functions. 


— t-err 


— p-err 
f 2 

— ^q-err 


flw- 


w-err 


(42) 
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Since the solution of equation (36) was based on small angle approximations and 
the angular motions required to maintain proper wirefeed may be quite large, the two sets 
of iterative equations (36) and (40) must be de-coupled by defining an artificial state 
variable, ©gr. Angle 0 ^’ is the angle required to maintain wirefeed control assuming the 
other constraints of position and downhand were met. The resulting iterative equation 
given in (43) and is evaluated subject to equation (44) being met. 


e j+1 =e j+ h* [f9 t _ err (9)J (43) 

[f9 t _ err (0) T ][K'][f9 t -err^)3 < epsilon (44) 
9 = [0i»02 / ©3 / 04 f ©5 » 06 / 07 f 08 / ®gl (45) 
© = [®2P 7 9 1P' 0 1R' 0 2R r 0 3R' 0 4R' 0 5R' 0 6R' 0 6R ' ^ < 46 ) 


The vector, 9, in equations (43) through (46) has nine components, 0 j through 9g, 
as defined previously and 09 is equal to the artificial state variable, 9^’. Matrix [J] is 
defined in equation (A-23). The vector function £9 t _ err is from equation (42), and finally, 
the matrix [K’j is a one dimensional extension of [K], the 8 x 8 matrix from equation (37). 
The added diagonal term K 99 controls the tolerance of the wirefeed orientation error. The 
results of a computer graphic simulation of this algorithm are now presented. 

GRAPHIC SIMULATION RESULTS 

The workcell layout, the six DOF robot, welding torch, the two DOF part positioner 
and the part were modeled using ROBOSIM [12], a NASA and Vanderbilt University 
developed computer graphic modeling package. Each mechanism is described using a 
procedure oriented instruction set. The completed model may be displayed on a variety of 
graphic displays including Tektronix 4014, a CAD system via IGES (Initial Graphics 
Exchange Standard), Evans & Sutherland PS330, GTI Poly 2000 or Silicon Graphics IRIS 


(other interfaces are planned). Figure ( 6 ) shows the workcell containing the robot, 
positioner and part. A close-up detail of the part is shown in Figure (7). The part for this 
example is a corrugated plate with the programmed weld seam as shown in Figure (7). 

Figures ( 8 a- 8 h) display the robot at several points along the path without downhand 
control in effect (K .77 and Kgg equal to zero in equation (37) with fixed positioner angles, 
0 jp and © 2 P)’ Th e assume d torch orientation is perpendicular to the surface tangent. 
Figures (9a-9h) show the same path points with downhand control (algorithm II) enabled. 
In Figures (9a-9h) we observe that the tangent to the weld path at the torch tip is 
horizontal, i.e., welding is occurring in the downhand position. An appreciation for the 
complex motion that must be performed by the torch can be gained by examining figure 
(10). Figure (10) depicts the torch trajectory with respect to the workcell reference frame 
during the downhand welding operation. This path bears little resemblance to the 
sinusoidal motion that would be observed in the part’s reference frame. The velocity profile 
of the torch tip also differs considerably when observed from the two reference frames. 
Figure (11) compares the torch’s velocity, observed in the workcell’s reference frame 
(dashed line), to the constant torch velocity observed in the part’s frame (solid line). We 
observe that the torch velocity must vary widely in workcell coordinates in order to 
maintain constant (. 10 ) units/sec seam velocity. 

Since algorithm III controls wire feed in addition to controlling the orientation of 
the part for downhand welding, a new part will be used to demonstrate this added feature. 
Figure (12) depicts this new part which is similar to the old part with a path that weaves in 
a sinusoidal manner as it travels over the surface. Algorithm III is demonstrated in Figures 
(13a-13h) which show an animation of the robot’s motion from the side at each time point. 
The front view clearly illustrates the control of downhand orientation while the 
corresponding top views in Figures (14a-14h) allow us to observe that the weld wire’s 
orientation is controlled. For clarity, the wire frame views in figures 14a through 14h have 
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been modified by rendering invisible, details that would obscure. Removed are the cell 
floor; positioner links zero and one; and robot links four and five. In the first frame of the 
animation sequence (14a) the weld wire orientation is not controlled since control is based 
on an approximation of the path tangent using a backward difference (equation 39). 
Subsequent frames, however, clearly show that the weld wire’s orientation is maintained. 

HARDWARE IMPLEMENTATION 

For those applications employing a CAD/CAM system to develop the path 
programs for the robotic arc welding cell, the transformation algorithm derived in this 
paper may be appropriately resident in the CAD/CAM system. It is only necessary for the 
weld assembly itself to be in the CAD/CAM display. The user specifies the position and 
orientation of the welding torch at each of the program points. The position and 
orientation is specified relative to the assembly. The kinematic transformation matrix, 
positioner axis commands, and velocity commands are downloaded from the CAD/CAM 
system to the robotic controller. 

For those applications employing a teach pendant to develop the path programs for 
the robot, the transformation algorithm would be resident in another computer. For each 
programmed point, the kinematic transformation matrix and the positioner axes values 
would be uploaded to this computer. The output of the transformation algorithm would 
then be downloaded to the robotic controller for welding. 

DISCUSSION AND CONCLUSIONS 

A transformation algorithm has been presented that permits the development of 
path programs for arc welding applications without the necessity of manually positioning 
the assembly to the downhand position or maintaining wire feed orientation for each point 
to be programmed. The algorithm automatically generates the required kinematic 
transformation algorithms, positioner axes command values and velocities to insure 
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downhand welding at constant weld travel speed with wire feed orientation controlled. The 
algorithm is compatible with CAD/CAM development of the path programs or with 
development of path programs via a teach pendant. As of this writing, the algorithm is 
being used to develop downhand adjusted path programs for a robot used in Space Shuttle 
Main Engine (SSME) weld process development [6-8] at the Marshall Space Flight Center. 
It is estimated that complex weld contours can be programmed in 75% less time by use of 
the transformation algorithm. On assemblies, such as the SSME, which contains over 2000 
welds. This represents a very substantial savings in costs and time. 

Since use of the transformation algorithm in developing path programs does not 
require that the weld assembly be positioned to the downhand position during the 
programming stage, it does not guarantee that the down position will be achievable for all 
welds programmed. It is hence important that this be established prior to welding. This may 
be done by means of a no-weld dryrun. However, a less time consuming method to check 
the path programs prior to welding is via computer simulation. 
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APPENDIX A: FORMATION OF THE JACOBIAN MATRIX 


Use of the Jacobian in manipulator control was first performed by Whitney [13] in 
his specification of resolved motion rate control (RMRC). RMRC has been applied to 
numerous man-in-the-loop manipulator control systems including the Space Shuttle’s 
Remote Manipulator System (RMS). The Jacobian is a matrix which relates joint angle 
rates to torch displacement and Euler rotation rates. The definition of the j-th column of 
the Jacobian matrix is given in equation (A-l) below. 


[J] 


element 


differential of 



torch-x in part frame 
w.r.t. the j-th joint 



torch-y in part frame 
w.r.t. the j-th joint 




torch-z in part frame 
w.r.t. the j-th joint 

(A-l) 

torch-x-rot. in part frame 
w.r.t. the j-th joint 



torch-y-rot. in part frame 
w.r.t. the j-th joint 



torch-z-rot. in part frame 
w.r.t. the j-th joint 


Joint angle definitions: 

for j=l, 8-s equals positioner angle second from base 
for j=2, ©:• equals positioner angle next to base 

for j=3, equals robot angle next to base 

for j=4, 8^ equals robot angle second from base 

for j=5, 8ij equals robot angle third from base 

for j=6 , 8j equals robot angle fourth from base 

for j=7, 0 j equals robot angle fifth from base 

for j=8, 0j equals robot angle sixth from base 


The physical interpretation of the Jacobian can be shown with the aid of figure (A- 
1). A wheel with axis of rotation w and a point (R) on its circumference with the vector r 
from the axis of rotation to the point (R). Both vector w and r are 3D vectors for this 
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example rather than the 4D homogeneous vectors used elsewhere in this paper. The 
instantaneous displacement of (R) due to a differential rotation of the wheel is given by the 
cross product of w and r. The instantaneous rotation of a reference frame located at (R) 
and fixed to the wheel would be w since the wheel is a single rigid body and as such all 
points in that body rotate about a fixed axis by definition. The Jacobian of this simple single 
degree-of-freedom mechanism is given in equation (A-2). 


[J] 


w x r 
w 


(A-2) 


The Jacobian can be similarly defined for the eight degree-of-freedom (DOF) 
mechanism formed by combining the two DOF positioner and the six DOF robot. 
Equations (A-3) contain this definition 


[J] 


— 1 —2 —3 —4 —5 —6 —7 —8 
w x w 2 w 3 w 4 w 5 w g w ? w 8 


(A-3) 


The end-effector displacements in part frame given by equations (A-4) to (A- 11): 

— 1 = -[Q]([G" 1 ][A 2p " 1 ](w X ([A lp - 1 ][Z p - 1 ]Z r T r (i)Er))) 

(A-4) 

v 2 = -[Q]([G _1 ][A 2p - 1 ][A lp - 1 ](w X ([Z p _1 ]Z r T r (i)Er) ) ) 


(A-5) 

v 3 = [Q]([G" 1 ][T p (i) _1 ][Z p - :L ]Z r (w X (T r (i)Er))) 

(A-6) 

Y 4 = [Q]([G" 1 ][T p (i)" 1 ][Z p _1 ]Z r A lr (w X (A 2r A 3r A 4r A 5r A 6r Er) ) ) 


(A-7) 
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~5 [ Q H[G ] [ Tp ( l ) 1 ]tZ p 1 ]Z r A lr A 2r (w X (A 3r A 4r A 5r A 6r Er) ) ) 



(A-8 ) 

~6 tQ](tG HTp(i) 1 ] [Z p 1 ]Z r A lr A 2r A 3r (w X (A 4r A 5r A 6r Er) ) ) 


(A-9 ) 

*7 - wms-'jtJpd)- 1 ] [^IViAAArt* x 

(AsrAgr^r) ) ) 


(A-10) 

~ I5]([G 1 ][Tp(i) 1 )[Zp 1 3 z r A lr A 2r A 3r A 4r A 5r(S 

X (Ag r Er) ) ) 

The end-effector rotations defined in part frame given by: 

(A-ll) 

*1 = -CQ3 C CG" 1 ] CA 2p “ 1 ]w) 

(A-12) 

-2 = ~[Q] ( [G" 1 ] [A 2 p _1 ] [Aj^p -1 ]^) 

(A-13) 

—3 = [Q3 ( [G -1 ] [T p (i) “ 1 ] tZ p _1 ]Z r w) 

(A-14) 

w 4 = [Q](tG- 1 ][Tp(i)-l][Z p - 1 ]Z r A lr ii) 

(A-15) 

*5 = tQHt G " 1 ][Tp(i)“ 1 ][Z p “ 1 ]Z r A lr A 2r w) 

(A-16) 

*6 = [QJCCG-^tTpfD-^fZp-ijz^^^^) 

(A-17) 

~ 7 ~ [Q H[G ^-JtTpfi) 1 ][Z p 1 ]Z r A lr A 2r A 3r A 4r w) 

(A— 18 ) 

~8 “ t Q HtG 1 ][T p (l) 1 ][Z p 1 ]Z^ ll ^ 2r A 3r A 4r A 5r w) 

(A-19) 

E = (0 0 0 1) T 
W = (0 0 1 0) T 

(A-20) 

(A-21) 

1 1 

o o o 
O O H 
OHO 
H O O 

II 

o? 

(A-22 ) 
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The Jacobian matrix defined in equation (A-3) has six rows and eight columns. The 
system of equations represented is under-determined and if equation (A-3) is to be used to 
control torch to part position and orientation, two additional constraints must be added. As 
stated in the text these additional constraints will require the part to be oriented so that 
welding occurs in the downhand position. These added constraints can be met by rotation 
of the part about a vector lying in the workcell’s horizontal plane. Equations (A-23) to (A- 
27) describe a Jacobian augmented with additional rotational vectors, dj and c^. These 
terms relate part rotations defined in workcell coordinates to joint angle rotations. The 
vectors dj and d 2 have been reduced to 2D vectors by the [R] operator which removed the 
x-component and the scale factor (fourth component). This was done since the x- 
component is vertical in the workcell reference frame definition in the example used. We 
further note that only columns one and two have non-zero d entries since only the two 
positioner angles are effective in orienting the part with respect to the gravity field. 


[J] 


1 

l< 

H 

—2 

—3 

>1 

in 

>1 

VO 

>1 

> 

>1 

1 

00 

>1 

— l 

—2 

*3 

—4 

—5 

*6 

—7 

—8 

d i 

—2 

0 

0 

0 

0 

0 

0 

— X 








Where d j and d 2 are given by: 


d l = [R][Z p ][A lp ]W 
d 2 = [R][Z p ]w 


[R] 

0 = 


0 10 0 
0 0 10 


[ 00 ] 


(A-23) 


(A-24 ) 
(A-25) 

(A-26) 

(A-27) 
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APPENDIX B: TORCH CONTROL WITH MINIMIZED JOINT DISPLACEMENTS 


The six degree-of-freedom position error function, f6p_ err , which was introduced in 
equation (15) vanishes as the error transformation of equation (12) becomes the identity 
matrix. This relationship is graphically depicted in figure (B-l). The ordinate of figure (B- 
1) is the norm squared of the function f£>p_ err corresponding to equation (21) while the 
abscissa is a one-dimensional representation of the function’s domain. The point, 0*, 
represents the solution of joint angles needed to maintain proper torch-to-part position. 
The point 0j represents the current estimate of the joint solution, 0 ; and 0j + j represents 
the next estimate of the joint solution. The derivation of an iterative equation using 
Newton’s Method begins with a truncated Taylor’s series approximation of the function, 
f6p. err as given in equation (B-l). 


— p-err + £S> 


— p-err ) + 


d— p-err 

de 


50 + H.O.T. 

(B-l) 


The vector, 50, represents the change in 0 between estimations, 0j and 0j + j, as 
given in equation (B-2). 


SO = ©j +1 - 0 j 


(B-2) 


Substituting (B-2) into (B-l) to solve for the 9j + 1 estimation we get equation (B-3). 


®j+l = + 


3f6_ 

— — p— err 

30 


-1 


— p-err (— + £6p_ err (0) 


(B-3) 


For our solution ft)p_ erT (fi + 56) is the zero vector, 0, and the definition of the 
Jacobian, [J], we may re-write this as equation (B-4). 
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©j+1 fij [J ] f6p_ err (0j) 


(B-4) 


] 

The matrix, [J], must be square to have an explicit inverse, equation (B-4) applies 
only for a six degree-of-freedom mechanism. In this case, however, the Jacobian has six 
rows and eight columns and cannot be inverted. Although [J] is non-square the relationship 
between [J] and variations in 0 and £6p_ err is still valid as expressed in equation (B-5). 

if6p_ er r = [J]i® (B-5) 

The system of equations defined in (B-5) is under-determined and an infinite 
number of solutions exist. One method [14] to overcome this problem is to impose the 
added constraint that the weighted sum-of-squares of joint angle displacements, 59, will be 
minimized. This constraint measure is given in equation (B-6). 

c = ^(£0 T W£0) (B-6) 

The matrix, W, is a diagonal weighting matrix to allow the the specification of 
factors which will increase the motion of preferred joint angles. A combined cost function 
is formed in (B-7) using a LaGrange multiplier, A, to append equations (B-5) and (B-6). 

s = ^ (ie T w6e) - [A T ][[J][ie] - f6 p _ err ] (B-7) 

The derivative of S with respect to £0 will be zero at the minimum condition 
yielding equation (B-8). 

ds m 

— = W50 - J L X = o (B-8) 

30 
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Solving (B-8) for 60 we get equation (B-9). 

69 = [W _1 ] [J T ]X (B-9) 

Multiplying (B-9) by [J] and substituting (B-5) into the result we get (B-10). 

[J]60 = [J] [W 1 ] [J T ]X = f6 p « err (B-10) 

Solving (B-10) for the LaGrange multiplier, X, and substituting into equation (B-8) 
we get equation (B-ll). 

60= [W _1 ][J T ] [[J][W _1 ][J T ]] [ fip-err 3 ( B “ 1:L ) 

Using our definition of 60 and 6f6 p err as given in figure (B-l) and equation (B-3) 
and (B-12) we get the resulting iterative equation in (B-13). 


££%>-err — " — p-err 

(B-12) 

e j+1 = 0 j - [w- 1 J T (JW- 1 J T )“ 1 ] f6 p _ err 

(B-13) 

— p-err T [ K 3 £^ p -err < epsilon 

(B-14) 

e j+1 = 0 j - J T (JJ T ) -1 ffi p _ err 

(B-15) 


Equation (B-13) is executed repeatedly until the weighted norm in equation (B-14) 
is less than the tolerance epsilon. For the case that equal weighing of joint displacement is 
desired joint angle weighting matrix, W, may be set to identity yielding equation (B-15). 
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APPENDIX C: COMPUTATION OF J' 1 

The eight by eight Jacobian matrix defined in Appendix A, equation (A-23) has the 
sparse form shown in equation (C-l). The inverse of the Jacobian is given in equation (C- 
2). We note that the inverse of the eight by eight matrix may be written in terms of the 
inverses of a six by six and a two by two matrix a result that significantly reduces the 
computation time. 


[J] 


A B 


C 0 


[j- 1 ] 


Where: 


0 

EC’ 1 ] 

[B -1 ] 

“[B _1 ]A[C _1 ] 


[A] = 



Hi —2 


[B] 


Y 3 v 4 v 5 v 6 v 7 v 8 
w 3 w 4 w 5 w 6 w ? w 8 


[C] = [ ^ d 2 ] 


(C-l) 


(C-2) 


(C-3) 


(C-4 ) 
(C-5) 


The matrix, B, in equations (C-l) and (C-2) is similar to the Jacobian of the six DOF 
manipulator alone, while the two by two sub-matrix, C, is a sub-matrix within the 
positioner’s Jacobian. The inverse of the eight by eight Jacobian (C-2) exists, if both the 
inverses of B and C exist. B'* exists when the robot’s Jacobian exists, and C‘* exists, for this 
example problem, when 9jp is not zero. 
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Figure 3. Robot path in part frame. 
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Figure 6. Simulated workcell. 






Figure 8. Robot path without downhand control. 
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Figure 11. Torch velocity. 
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Six DOF error function characteristic. 
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